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ABSTRACT 

The structure of the hot downstream region below a radiative accretion shock, such as 
that of an accreting compact object, may oscillate due to a global thermal instability. 
The oscillatory behaviour depends on the functional forms of the cooling processes, the 
' energy exchanges of electrons and ions in the shock-heated matter, and the bound- 

ary conditions. We analyse the stability of a shock with unequal electron and ion 
temperatures, where the cooling consists of thermal bremsstrahlung radiation which 
promotes instability, plus a competing process which tends to stabilize the shock. The 
effect of transverse perturbations is considered also. As an illustration, we study the 
special case in which the stabilizing cooling process is of order 3/20 in density and 
5/2 in temperature, which is an approximation for the effects of cyclotron cooling in 
. magnetic cataclysmic variables. We vary the efficiency of the second cooling process, 

' the strength of the electron-ion exchange and the ratio of electron and ion pressures 

Q> , at the shock, to examine particular effects on the stability properties and frequencies 

0^ ' of oscillation modes. 

' Key words: accretion — shock waves — stars: binaries: close — stars: white 

dwarfs 

o 

C3 1 INTRODUCTION 

^ ■ Depending on the form of cooling processes present, radiative shocks may suffer thermal instabilities. The resulting sequence 
' of oscillatory modes has frequencies and stability properties determined by the functional form of the energy loss terms, 
. the rate of energy exchange between ions and electrons, and the particular boundary conditions of the system. Oscillations 
' modulating the structure of the hot downstream column cause variations in the observable emissions. Radiative shocks in 
several different astrophysical systems show thermal instability, for example the interactions between supernova shocks and 
the interstellar medium, and accretion streams onto magnetic and non-magnetic compact objects. Many investigations have 
considered the stability properties of radiative shocks, either through numerical calculations or perturbative analyses (e.g. 
Falle 1975, 1981; Langer, Chanmugam & Shaviv 1981, 1982; Chevalier & Imamura 1982; Imamura, Wolff & Durisen 1984; 
Chanmugam, Langer & Shaviv 1985; Imamura 1985; Bertschinger 1986; Innes, Giddings & Falle 1987a, b; Gaetz, Edgar & 
Chevalier 1988; Wolff, Gardner & Wood 1989; Imamura & Wolff 1990; Houck & Chevalier 1992; Wu, Chanmugam & Shaviv 
1992; Toth & Draine 1993; Dgani & Soker 1994; Imamura et al. 1996; Saxton, Wu & Pongracic 1997; Saxton et al. 1998). 

In accreting white dwarfs the flow meets a stand-off accretion shock where the inwards falling matter is suddenly decel- 
erated to subsonic settling flow. The shock sits above the white dwarf surface at a height determined by the path length of 
the radiative cooling of the downstream material a^s jvgtcooi, which depends on the free-fall velocity of the white dwarf 
Wff = (2GMwd/J?wd)^^^, and on the cooling timescale tcooi ~ k^T/A, where A is a radiative cooling loss function. In cataclysmic 
variables (see Warner 1995), the white dwarf accretes mass from a close red dwarf companion star. Thermal bremsstrahlung 
cooling is an important cooling process in the post-shock region. In cases where the magnetic field of the white dwarf is strong, 
(Bwd ^ 10 MG), (e.g. Lamb & Masters 1979; King & Lasota 1979), cyclotron radiation provides another important cooling 
process. These competing processes have different dependences on the local flow variables, and therefore they influence the 
shock stability and oscillatory behaviour differently. 



Langer et al. (1981) carried out time-dependent numerical investigations of accretion onto a white dwarf in a planar 



geometry. Their results show quasi-periodic oscillations due to thermal ins tabili ty when the cooling is via bremsstrahlung 



radiation. In the linear perturbative analysis as in Chevalier & Imamura (1982), the radiative shock with bremsstrahlung 



cooling was shown to have a stable fundamental mode and unstable overtones. They also investigated other cases with 
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a power-law cooling function in density p and temperature T, A oc p"' ^^ with various choices of the power indices (e.g. 
(a, b) = (0.5, 2) for bremsstrahlung cooling). These investigations show that radiative shocks where the index b is higher are 
more strongly sta bilized against perturbations. 

Bertschinger (198(: ) studied the effect of placing the radiative shock in a spherical geometry. Again the cooling function was 
a single power-law of density and temperature. In addition to radial perturbations, transverse perturbations were investigated, 
which are expressed in terms of a scaled transverse wavenumber in addition to the usual oscillation frequency. The modes 
that are stable in the purely radial analysis of Chevalier & Imamura (1982) are destabilized in the presence of transverse 
oscillation, but all modes are stable in the limit of large wavenumber. 

The stability of physically extended shocks was studied by Houck & Chevalier (1992) for radial accretion onto compact 
objects, with a Newtonian gravitational potential. The adiabatic index of the gas and the indices of the single power-law 
cooling term were varied. In the limit of a non-extended shock envelope their calculations yield the same results as Chevalier 
& Imamura (1982). However even a slight spherical extension of the post-shock region has an effect on the stability: stabilizing 
the first and second harmonics but destabilizing the higher order modes. 

Physical factors other than the cooling processes and geometry also affect the stability properties of a radiati ve sh ock. 
The effects of mass loss from the sides of the post-shock accretion column were examined by Dgani and Sokcr (1994) by 
including a sink term in the mass continuity equation. They found that the transverse leakage has a stabilizing effect in the 
presence of several single power-law radia tive c ooling processes, and this stabilization was less effective when the cooling had a 
lower temperature index. Toth & Draine (1993) investigated radiative shocks that were subject to thermal instability but were 
supported by a transverse magnetic field. This magnetic pressure support is decisive in determining the stability properties, 
rather than the characteristics of the radiative cooling alone. 

The consequences of inequality between electron and i on tem peratures, and the effect of the energy exchange between 
these fluid components, was considered by Imamura et al. ( [1996| ), with the effec ts of bremsstrahlung and Compton cooling. 
Their perturbations included a transverse component, as in Bertschinger (1986), although the geometry of the system was 
planar. Modes with transverse components were less stable than purely longitudinal equivalents, and the presence of the 
electron-ion exchange also promoted instability. These new effects were qualitatively similar under different choices of the 
temperature index of the single cooling process. 

In these studies, the usual treatment of the radiative cooling processes was to use an approximate energy loss term that is 
a single power-law in local fluid properties. In the presence of multiple cooling processes, intermediate power-law indices tend 
to be used. This is not an adequate treatment for the detailed effects of the competition between cooling processes, which is 
especially important when one process, such as thermal bremsstrahlung radiation, promotes thermal instability, while another 
process, such as cyclotron emission, has a stabilizing influence. The complicated competition between such processes must be 
treated by taking an ex plicit sum of independent power-law contributions for the respective processes. 

Chanmugam et al. (1985) performed numerical calculations of the accretion onto a magnetic white dwarf with a sum of 
separate bremsstrahlung and cyclotron power-law contributions. When cyclotron cooling was strong, the shock was stabilized 
against oscillations. Further examinations (Wu et al. 1996) revealed that the stabilizing influence depends on the magnetic field 
strength. In a weak-field regime {B ^ lOMG) the competition of the bremsstrahlung and cyclotron cooling, each dominating 
in different phases of a two- phase oscillatory cycle, allows the perpetuation of small-amplitude oscillations (Wu et al. 1992). 

Hujeirat & Papaloizou (1998) performed one- and two-dimensional numerical calculations to investigate the structure and 
time- dependent behaviour of accretion onto compact objects in the presence of bremsstrahlung cooling, and with consistent 
solution of the radiative transport and MHD equations. Magnetic field strengths were investigated up to ~ lO"'^ times the 
values found in magnetic cataclysmic variables (mCVs). They found that radiation- matter coupling in the flow below the 
photosphere reduces the amplitudes of the shock oscillations, but that the oscillations are sensitive to the lower boundary 
condition. In weak-field cases, latitudinal flows out of the accretion column reduce or suppress oscillations; and in stronger 
field cases there remains some transverse kink s in the magnetic field near the shock, which tends to suppress oscillations 
because of the effect found by Toth & Draine (1993). 

A composite coolin g function was used in the single-temperature perturbative stability analyses of Saxton et al. (1997) 
and Saxton et al. (1996 ). The energy loss term consisted of a sum of terms for thermal bremsstrahlung radiation (Abr oc p^T"'^) 
and a second power- law process with a destabilizing influence (A2 oc pT*" for 6 <: 1). The cases considered included that of 
an effective cyclotron cooling term (Acy oc p0-i5y2.5-j appropriate for the conditions and flow geometry of accreting magnetic 
white dwarfs. The analyses examined effects of different values of the efficiency of the second cooling process, which depends 
on the magnetic field strength in the cyclotron-cooling case. Increasing the cyclotron cooling strength stabilized each mode 
monotonically, but this happened in a more fundamentally complicated manner than straightforward comparison of cooling 
and oscillatory timescales would suggest. The order and manner in which the harmonics stabilized depends on the indices of 
the second cooling process. 

In the paper we generalize the previous work in Saxton et al. (1997) and Saxton et al. (1998) for bremsstrahlung-cyclotron 
radiative accretion shocks in mCVs, by including the effects of unequal electron and ion temperatures and the possibility of 
transverse perturbations (for a corrugated shock). The equations for a radiative post-shock accretion flow are presented in 
section ^ and the terms and formulation of the perturbation analysis are developed in section ^. In section ^ we calculate 
and interpret the frequency and stability properties of the eigenmodes, and explore the effects of varying the parameters of 
the two-temperature condition, and the relative strengths of the competing cooling processes. In section ^ we discuss the 
significance of our results in contrast with previous studies dealing with cases with different physical processes or which were 
less general. We conclude in section ^. 
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2 FORMULATION 



We assume the adiabatic index 7 = 5/3 for an ideal gas, and the equation of state 

P _ pKT 

where fcg is the Boltzmann constant, and is the mass of the a hydrogen atom. The time-dependent mass continuity, 
momentum and energy equations for the planar post-shock accretion flow are: 



(1) 



+v - v + VP = 0, 
(|+v.v)p-7^(|+v.v)p = -(7-l)A, 

+ v-v)p=(7-l)(r-A), 



P / d 
P 

Pc / d 
7- 



(2) 
(3) 
(4) 
(5) 



where P, Pe, v and p are respectively the total pressure, electron pressure, velocity and density of the stream; A is the 
composite radiative cooling function, and T is the electron-ion energy exchange. 



4v27re rieTii InC 



(me/mi)9e 



(6) 



where InC is the Coulomb logarithm, c is speed of light, e is the electron charge, 9^ = kj^Te/nieC^ and 9i = k^Ti/m\c^ are 
the electron and ion temperatures (see, e.g. Melrose 1986). The electron and ion number densities are rie = Zp/{Zme + mi) 
and rii = p/{Zmc + mi) if the ion charge is Z (Z — 1 for completely ionized hydrogen plasma). 

Optically thin thermal bremsstrahlung radiation provides one of the two cooling processes present, and the other is 
taken to be a single power-law cooling term. To simplify the analysis, we express the total cooling function in terms of the 
bremsstrahlung cooling term and a function expressing its functional form compared with the bremsstrahlung cooling term. 
The parameter eg is the efficiency of the second cooling process compared to bremsstrahlung cooling, evaluated at the shock, 
which depends non-linoarly on the white dwarf's magnetic field strength, mass and radius, and the pre-shock density and 
cross-sectional circumference and area of the accretion flow. 



A^ 



Kr + A2 = Abr 



1-hes 



Pc,. 



~I3 



(7) 



where ps and Pe,s are the density and electron pressure at the shock, Abr = Cp {Pe/p) and the constant C = 
{2T^k^/ZmeY^'^{2^TTe^ /?>hmeC?){ix/k^'m^Y^^g^, with h being the Planck constant, mc the electron mass, mp the proton 
mass, p. the mean molecular weight of the gas and the Gaunt factor (see Rybicki & Lightman 1979). The numerical value 
of C is 3.9 X 10^® in c.g.s. units, assuming that p = 0.5 and « 1. The indices a and (3 are constants describing the functional 
form of the second cooling process; they adopt particular values for radiative accretion shocks in different physical contexts, 
but we retain them as parameters for the sake of the general formulation, until results for mCVs are outlined below. (In 
general for a cooling function A2 oc p°'T^, a = 6 — 1/2 and /? = 3/2 — a -|- 6.) 

Expressing the composite cooling function in terms of the bremsstrahlung cooling assures us that the second cooling 
function goes to zero as the density becomes infinite. Using this form assures us that the stability results depend only on 
the indices of the cooling function and are insensitive to the implementation of the boundary condition at the white dwarf 
surface, which becomes unambiguous. Alternatively, the cooling could be stopped at a finite density or temperature, with a 
softened form of the bremsstrahlung term and a different choice of lower boundary. However this approach is mathematically 
equivalent to the one we have taken except that the new boundary conditions are less obvious, especially for the perturbed 
variables. 



3 PERTURBATION ANALYSIS 

The shock position Xg and velocity Vs are subject to a first-order perturbation (as in Imamura et al. 1996): 

Vs = Vsie (8) 

Xs = a::sO + x^ie^^^'^* . (9) 

with frequency (jJ and a transverse wavenumber k allowing a transverse component for corrugated oscillation. Subscripts and 
1 correspond to the stationary state solution and perturbed quantities respectively. In the stationary solution, the shock is at 
rest: Vao = 0. The longitudinal perturbation of the shock position and velocity are related by Va\ = XsiuJ- The frequency of 
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the perturbation is a;, and k is the wavenumber of a transverse component of the perturbation in the y direction. Eigenmodes 
are labelled by a transverse wavenumber and dimensionless eigenfrequency, which are defined as: 



= U) 



(10) 
(11) 



The eigenfrequency is complex, 5 = 5^ + iS-^, where 5j is the oscillatory part proportional to the frequency of the particular 
mode. Stability is indicated by the sign of the growth/decay term (5j^. Positive values indicate instability; negative values 
indicate stability. The physical scales of the system are vs the free-fall velocity of the white dwarf, pa the accretion density, 
and Xso the equilibrium shock height. The size of the perturbation is indicated by the parameter e = Vsi/vff, and therefore 
Xs\ — exao/5. 

To simplify the calculations, the hydrodynamic quantities are expressed in terms of dimensionless variables, allowing the 
scales of the system to be removed from later equations. These dimensionless variables are expressed in terms of the stationary 
solutions and perturbations of the velocity, density, total pressure and electron pressure, varying with position in the accretion 
column ^ = x/xs- The white dwarf surface is at ^ = and the shock is ^ = 1. 

Pa-Co(C)(l+eAc(C)e'^^+"') (12) 

-va ■ (ro(C) (1 +£A.(e)e''=''+-*) , ero(OA,(C)e"'=^+"') (13) 

Pa4-vro(C)(l + eA4C)e"'"+"') (14) 

P.vl ■ 7reo(C) (l + eAe(C)e"=^+"') . (15) 

where the complex functions A;;, At, \y, A^ and Ac represent the effect of the perturbation on the downstream profiles of 
density, vertical velocity, transverse velocity, total pressure and electron pressure respectively. 

The space and time derivatives of a flow variable /, in a frame following the shock, are given by: 

iky + ujt ^ 



£1 

Dx 

£1 
Dy 

£1 

Dt 



1 

XbO 

21. 

dy 

21 

dt 



1 

a;so 

Vft 
XbO 



■ £.1^ iky+ujt 

le—e 





21 
21 
21 



(16) 
(17) 
(18) 



None of the hydrodynamic variables has an explicit dependence on position within the accretion column. These shock-frame 
derivatives are used to expand the equations for mass, momentum and energy continuity, whence are extracted expressions 
describing the stationary solution and a set of differential equations for the perturbed variables. 

We consider expressions in terms of the dimensionless stationary-solution velocity (ro = —vq/vs) to simplify the form of 
the equations. The stationary solution is obtained by solving the hydrodynamic equations without the time-dependent terms, 
and it is completely described by two simple algebraic equations for the density and pressure, 



Co = 1/ro 

TTO = 1 — To 

plus two differential equations: one for the velocity-position proflle and the other for the electron pressure: 
dj _ 77ro - tq 



d-Ka 1 

dro ro 
where 



1 — — (77ro — ro) — 77rc 
A 



(7- l)(Pai'fi/a:so) 



(7 - l)'/'cV'c 



■ ro ■ 



27re 



and 



A = (7 - l)(paWff/Seo)"'A = (7 - l)^/>c 



[l + es/(ro,^c)l 



(19) 
(20) 

(21) 
(22) 

(23) 
(24) 



are appropriate dimensionless forms of the electron-ion exchange and total cooling function respectively. The parameter i/)ci 
describes the speed of the electron-ion exchange compared to the cooling; when ijjd ^ 1 the two-t empe rature treatment is 
justified; otherwise the system behaves like the single-temperature shock described in Saxton et al. ( 1995 ) . In physical terms 
the constants tpc and tpd are: 

XsOpa. 



■c, 



(25) 
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and 



(26) 



and T(crs) ~ 1 (see Saxton 1999). The function /(ro,7re) describes the local ratio of the power of the second cooling process 
relative to the bremsstrahlung emission. It takes the form: 



/l+Os\" a 



(27) 



Pc/Pi)B is the ratio of electron to ion pressures at the shock surface. Upon substitution of these expressions, (|l9l 



and U a = (f „ , . , „ 

|2^ , p^ , p7| ) the equations for the stationary solution velocity (^) and pressure profiles ( |22| ) become: 
7 — (7 + l)ro 



dro (7 - l)7/)c [1 + es/(-ro,7rc)] 



dro 



1 

To 



1 - To - 27rc 



1 + es/(ro, TTc) 



roTTc 



(7 — (7 + 1) To) — 77rc 



(28) 
(29) 



For the case of bremsstrahlung cooling alone (es = 0) there exists an analytic solution (Aizu 1973). For the general two-process 
cooling function it is necessary to carry out a numerical integration. 

The hydrodynamic equations expressed in terms of the dimensionless perturbed variables give rise to a set of five coupled 
complex first-order differential equations, which can be expressed in matrix form as: 



(30) 



where a prime denotes derivatives in ^ and the functions on the right hand side collect the terms which lack derivatives of the 
perturbed variables. In general, for a cooling function composed of a sum of several power-law terms, they are: 



-To 


— To 








" 




r K 1 









— To 





TO - 1 







K 




F2 








-To 














F3 


7To7ro 








— TOTTO 







\„ 




Fi 


7T07re 











— ToTTo 








F5 



Fi{to,7Vc,0 = 

F2{ro,TVc,0 = 

F3{ro,TVc,0 = 

F4,iro,-rvc,0 = 

F5{To,nc,0 = 



-^(Inro)' — 5A<; + iKToXy 

-{5 — tq)Xt + ^(Inro)' + ToA(; — tq{X^ — Ar) 

-(5 — To)Ay + i/t(l — To)A,r + inTQ^^/S 

-TTnS{\^ - 7A^) 



A 



3 1 1 ^ 

Xff2(To,7rc)A<; + -gi{To,TVc)Xc - Xt - X^ + 

2 2 To 



— 7re<5(Ae — 7A1;) — A 
+ t 



3 1 1 ^ 

--52(T0,7rc)Ac + -3l(T0,7rc)Ac - Ar - A^ + - - 



To 



5 3 TToA^ 

2 ^ 2 TTO 



- 2^cAe 1 C 

""0 At — Ac -I- J 

27rc To 



where the functions gi{To,Trc) and (72(To,7rc) are defined as: 
2esa/(ro,7rc) 



5l(T0,7I"e 



l-f 



1 -I- es/(T0, TTc 



ff2(T0,7rc) = 1-3 



es/3/(To,7rc) 



1 + es/(T0, TTc 



(31) 
(32) 
(33) 

(34) 



(35) 

(36) 
(37) 



The perturbed transverse velocity variable Aj, has an inconspicuous infiuence. If not for its presence in the Fi function, 
Xy would be decoupled from the other perturbed variables, and could be integrated separately after finding the solutions in 
the other variables. 

The matrix in (^^ is non-singular everywhere in the post-shock region except at tq = 0. In a practice this singularity 
doesn't matter because the numerical integration can be taken to some small value ro = A > and the results reach a limit 
as A 0. The matrix equation is inverted using (Ell) to yield equations for the perturbed variables. 



d 

dTO 







1 


-1 





1/to 







" Fi ' 


Xr 


1 

^ I 


-77ro/ro 


1 





-1/to 







F2 


Xy 








-(77="o - To) /to 










F, 


Att 


7 


-7 





I/tto 







Fi 


Ac . 




7 


-7 





7/to 


-(77^0 - To)/To7rc 




Fs 



(38) 



These resulting differential equations are split into real and imaginary parts to yield ten first-order differential equations in 
real eigenfunctions, A^^, A^^, Xr-^, At,, ^y-^y -^Hji -^ttr, A^,, Ae^^, and Ac,, plus the known parameters and variables of the 
static solution, q, /3, ^(to) and 7rc(To). This decomposition is straightforward because only the vector parts of (Bq) are complex. 
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To obtain the equations for real variables, we substitute the real parts of the A-variables in the left side and real parts of the 
F-functions in the right side. Similarly, the equations for the eigenfunctions' imaginary parts are obtained by substituting 
imaginary parts of the A-variables and F-functions. 

Optically thick cyclotron emission provides an important cooling process for accretion onto strongly magnetic white 
dwarfs. An optically thick radiative process cannot in general be described by a simple local energy loss function. A full 
treatment requires simultaneous solution to the time-dependent hydrodynamic and radiative transfer equations. However 
for the particular geometry of the field-aligned accretion flow in mCVs and the physical conditions typical of such systems, 
cyclotron radiation is optically thick up to a critical frequency at which it becomes optically thin. The total power escaping 
from the post-shock region locally in the accretion column approximately takes the form of a power-law cooling term Acyc oc 
^3/20ji^5/2^ which in the terms of corresponds to indices a = 2.0 and /3 — 3.85 (see Langer et al. 1982, Wu, Chanmugam 
& Shaviv 1994, Cropper et al. 1999). This approximation renders the linear analysis tractable. 

For the special case of a mCV radiative shock with bremsstrahlung and cyclotron cooling present, the parameter eg is 
related to the magnetic field at the accretion pole. The expression for eg is as described in Wu (1994). In the single-temperature 
limit (Ts = 1, but the value in the two-temperature conditions depends on the conduction of energy upstream of the ion shock, 
carried by the electrons. Our analysis treats CTs as a parameter. 



4 RESULTS 

4.1 Eigen-modes 

The ten real differential equations are integrated numerically using a Runge-Kutta method. In terms of the dimensionless 
perturbed variables, the boundary conditions at the shock are tq = 1/4, ^ = 1, ttc = |(1 -I- 1/(Ts)^^, \q — Q, \r — —3, 
\y = Sin/S and A,r = Ao = 2. At the dwarf surface (to = 0, ^ = 0, ttc = 1/2) there are no specific conditions on the values 
of the perturbed density and pressures, but the stationary wall condition At = applies (see Chevalier & Imamura 1982; 
Imamura et al. 1996; Saxton et al. 1997). Integration proceeds between ro= and ro = l/4 for trial values of 5^ and 5^. Values 
of the S's are sought which yield A^ = when the differential equations are integrated to tq = (using the search method 
described in the appendices of Saxton et al. 1998) Each of these eigenvalues corresponds to an oscillatory mode of the shock. 
The modes form an indefinite sequence consisting of a fundamental mode and a succession of overtones. 

The survey of the complex frequency eigenplane is conducted in a region of S appropriate for the lowest modes. The range 
of frequency explored is from to 3.5; the growth/decay term is taken from 0.4 down to —0.4 or lower if necessary in 
cases of higher k. The perturbed velocity eigenfunction is evaluated at the white dwarf surface for each sampled of 5. Values 
of — log jA^lwd are plotted and the maxima of this quantity indicate the 5 eigenvalues. 

This procedure is performed for single-temperature (cts = 1, large i/)ei) and two-temperature (cts < 1, tpci finite) cases. 
Eigenfrequencies of the six lowest modes are calculated under the conditions as — 0.2,0.5, 1 and -i/joi = 0.1,0.5, 1 to test the 
effect of the unequal electron and ion pressures at the shock, and the rate of electron-ion energy exchange in the accretion 
column. For each case we investigate the modes' frequencies and stability properties under conditions with bremsstrahlung 
cooling alone (ts = 0), cyclotron cooling comparable to the bremsstrahlung (es — 1), and cyclotron cooling dominant (eg = 100), 
as evaluated at the shock surface. These eigenvalues are tabulated in Table |l]. 

Applying the special restricted choice of the eg = exactly recovers earlier results of Imamura et al. (1996) for a two- 
temperature shock with only a single bremsstrahlung cooling process. Taking the limits of high tpd and as — 1 reduces our 
equations to reproduce the single-temperature case with multiple cooling, as in Saxton et al. (1997). Taking this limit and 
Ea = together reproduces the single-temperature bremsstrahlung-cooling only results of Chevalier & Imamura (1982). 

The frequencies and stability behaviour of the fundamental and first overtone in the presence of the two-temperature 
conditions and composite cooling was also investigated for non-longitudinal perturbations. This corrugation of the shock is 
represented by non-zero values of the dimensionless transverse wavenumber k. 

4.2 Frequency quantization 

Under the simple conditions, for either bremsstrahlung-dominated cooling or a single-temperature shock, the oscillatory part 
of the eigenfrequency follows a quantized sequence like the modes of a pipe open at one end. In a linear fit to this sequence, 
the harmonics n follow ~ 5jQ(n— 1/2) + 5^, with frequency spacing 5^^ and a small constant correction offset S^-.. The two- 
temperature parameters as and tpai and the secondary process' cooling efficiency es determine the values of these constants. 
Parameter fits were performed for the first six harmonics under different conditions, and the results for the constant and 
linear coefficients are shown in Table ^. 

For a purely bremsstrahlung-dominated shock the mode spacing (5jq ~ 0.6 regardless of the two-temperature conditions. 
As es is increased under given as and conditions, 5^^ becomes steadily smaller. For very efficient cyclotron cooling (es = 100) 
the spacing is reduced to 5^^ ~ 0.3. The decline of Sj^ with increasing es is fastest for cases with high ipd; and increasing 
ipei for given {as,es) (approaching single-temperature conditions) makes 5jq smaller. In other words, when the electron-ion 
exchange is weak, the effect of es on the mode frequencies is reduced. 

The frequency offset constant 5^ increases steadily with increasing es for all observed values of the two-temperature and 
cooling parameters ((Js,i/'ci). For an increasingly cyclotron-cooling dominated flow in the presence of strong two-temperature 
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Table 1. Eigenvalues 5^, 5^ in the absence of transverse perturbation (k = 0), for various conditions of electron/ion pressure ratio at 
the shock {<7b), strength of electron-ion energy exchange (V'ei)) ^ind relative efficiency of cyclotron cooling (es). 





V'ei 




= 




= 1 


Es = 


100 


0.2 


0.1 


-0.001 


0.307 


-0.057 


0.343 


-0.113 


0.335 






0.048 


0.902 


-0.040 


0.808 


-0.073 


0.567 






U.UDO 


1 CCOQ 

l.OZif 


— U.UoU 


l.oOo 
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2.444 


-0.036 


1.758 






0.114 


3.385 


0.027 


2.995 


-0.046 


2.154 


0.2 


0.5 


-0.007 


0.306 


-0.106 


0.340 


-0.166 


0.247 






0.047 
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-0.043 


0.528 
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0.085 
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-0.040 


1.756 


-0.088 


1.104 






0.087 


2.739 


-0.054 


2.253 


-0.071 


1.406 






0.106 


3.349 


-0.026 


2.770 


-0.086 


1.710 


0.2 


1.0 


-0.009 


0.306 


-0.123 


0.332 


-0.163 


0.211 






0.047 


0.891 


-0.056 


0.744 


-0.056 


0.520 
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n '7ni 

U. (Ul 






0.085 


2.115 


-0.046 


1.717 


-0.109 


1.006 






0.086 


2.733 


-0.067 


2.193 


-0.081 


1.271 






0.105 


3.342 


-0.043 


2.701 


-0.094 


1.541 


0.5 


0.1 


0.003 


0.304 


-0.025 


0.322 


-0.057 


0.309 






0.051 


0.898 


0.006 


0.849 


-0.027 


0.582 
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U.U14: 
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2.132 


0.044 


2.002 


0.003 


1.416 






0.096 


2.756 


0.044 


2.585 


0.002 


1.819 






0.114 


3.370 


0.064 


3.161 


0.008 


2.232 


0.5 


0.5 


-0.006 


0.306 


-0.065 


0.332 


-0.135 


0.264 






0.048 


0.892 


-0.019 


0.823 


-0.023 


0.551 






U.UOi 
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— U.UID 


i.ooy 


— U.Uol 


U.C5DO 






0.086 


2.117 


0.011 


1.930 


-0.049 


1.218 






0.088 


2.735 


0.006 


2.490 


-0.051 


1.547 






0.107 


3.345 


0.031 


3.047 


-0.059 


1.906 


0.5 


1.0 


-0.008 


0.305 


-0.078 


0.332 


-0.148 


0.239 






0.048 


0.891 


-0.023 


0.814 


-0.036 


0.549 






U.UOl 


l.OUo 


— U.Uzo 


l.O ( O 


— u.uyz 


U.ttlD 






0.085 


2.113 


0.006 


1.910 


-0.073 


1.159 






0.087 


2.730 


-0.002 


2.463 


-0.064 


1.448 






0.106 


3.339 


0.023 


3.015 


-0.083 


1.781 


1.0 


0.1 


0.010 


0.299 


-0.001 


0.303 


0.002 


0.277 






0.056 


0.892 


0.036 


0.856 


0.027 


0.564 
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U.U4 1 
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0.093 


2.118 


0.064 


2.025 


0.052 


1.378 






0.095 


2.735 


0.064 


2.610 


0.043 


1.798 






0.112 


3.345 


0.078 


3.192 


0.060 


2.199 


1.0 


0.5 


-0.005 


0.304 


-0.039 


0.319 


-0.096 


0.262 






0.049 


0.890 


0.009 


0.846 


0.007 


0.554 






0.063 


1.507 


0.019 


1.430 


-0.032 


0.909 






0.087 


2.112 


0.042 


1.995 


-0.000 


1.257 






0.089 


2.729 


0.040 


2.576 


-0.022 


1.610 






0.107 


3.337 


0.061 


3.150 


-0.014 


1.992 


1.0 


1.0 


-0.007 


0.305 


-0.050 


0.322 


-0.122 


0.247 






0.049 


0.890 


0.004 


0.843 


-0.015 


0.559 






0.062 


1.506 


0.012 


1.425 


-0.059 


0.876 






0.086 


2.110 


0.038 


1.988 


-0.025 


1.230 






0.088 


2.727 


0.035 


2.566 


-0.038 


1.545 






0.106 


3.335 


0.056 


3.139 


-0.046 


1.915 
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Table 2. Mode frequency spacing 5jq and offset 5^ for a linear fit 5, (5jq (n — 1/2) + &^ , under different two-temperature conditions 
(cTsiV'ei) S'lid efficiency of the second cooling process (eg). 



0"s 








Sc 








^lO 


<5c 


0"s 






5io 


5c 


0.2 


0.1 





0.617 


-0.013 


0.5 


0.1 





0.615 


-0.014 


1.0 


0.1 





0.610 


-0.015 


0.2 


0.1 


1 


0.534 


0.039 


0.5 


0.1 


1 


0.571 


0.014 


1.0 


0.1 


1 


0.580 


0.000 


0.2 


0.1 


100 


0.374 


0.076 


0.5 


0.1 


100 


0.392 


0.051 


1.0 


0.1 


100 


0.391 


0.028 


0.2 


0.5 





0.610 


-0.010 


0.5 


0.5 





0.609 


-0.011 


1.0 


0.5 





0.608 


-0.011 


0.2 


0.5 


1 


0.489 


0.057 


0.5 


0.5 


1 


0.547 


0.030 


1.0 


0.5 


1 


0.569 


0.013 


0.2 


0.5 


100 


0.294 


0.080 


0.5 


0.5 


100 


0.330 


0.069 


1.0 


0.5 


100 


0.348 


0.055 


0.2 


1.0 





0.609 


-0.011 


0.5 


1.0 





0.608 


-0.010 


1.0 


1.0 





0.608 


-0.011 


0.2 


1.0 


1 


0.476 


0.060 


0.5 


1.0 


1 


0.540 


0.032 


1.0 


1.0 


1 


0.566 


0.015 


0.2 


1.0 


100 


0.263 


0.086 


0.5 


1.0 


100 


0.307 


0.077 


1.0 


1.0 


100 


0.333 


0.063 



Table 3. Parameters of the frequency spacing for a quadratic fit, 5j Ri ^q'"^ + 5jo(" ~ 1/2) + under different two-temperature 
conditions {o-s,i/!ei) and efficiency of the second cooling process (eg). 









5io 










ts 


5io 


5c 




Os 




Cs 


5io 


5c 




0.2 


0.1 





0.604 


-0.002 


0.002 


0.5 


0.1 





0.602 


-0.033 


0.002 


1.0 


0.1 





0.601 


0.007 


0.001 


0.2 


0.1 


1 


0.484 


0.080 


0.007 


0.5 


0.1 


1 


0.541 


0.038 


0.004 


1.0 


0.1 


1 


0.564 


0.013 


0.002 


0.2 


0.1 


100 


0.291 


0.144 


0.012 


0.5 


0.1 


100 


0.319 


0.113 


0.011 


1.0 


0.1 


100 


0.327 


0.082 


0.009 


0.2 


0.5 





0.596 


0.001 


0.002 


0.5 


0.5 





0.595 


0.001 


0.002 


1.0 


0.5 





0.595 


0.000 


0.002 


0.2 


0.5 


1 


0.436 


0.102 


0.007 


0.5 


0.5 


1 


0.508 


0.062 


0.005 


1.0 


0.5 


1 


0.541 


0.036 


0.004 


0.2 


0.5 


100 


0.252 


0.115 


0.006 


0.5 


0.5 


100 


0.279 


0.112 


0.007 


1.0 


0.5 


100 


0.292 


0.101 


0.008 


0.2 


1.0 





0.594 


0.002 


0.002 


0.5 


1.0 





0.594 


0.001 


0.002 


1.0 


1.0 





0.593 


0.002 


0.002 


0.2 


1.0 


1 


0.428 


0.100 


0.007 


0.5 


1.0 


1 


0.501 


0.065 


0.006 


1.0 


1.0 


1 


0.536 


0.041 


0.004 


0.2 


1.0 


100 


0.245 


0.100 


0.003 


0.5 


1.0 


100 


0.282 


0.098 


0.004 


1.0 


1.0 


100 


0.298 


0.093 


0.005 



eflFects, 5^ — > ^(5jo' which means that the nature of the modes becomes more analogous to a pipe with two open ends. The 
increased importance of the two temperature effects in the large Es regime means that the fixed wall boundary condition loses 
importance and the global stability of the accretion column depends more on the disparity of electron and ion temperatures 
within the stream than on the radiative cooling processes (cf. Imamura et al. 1996). The frequency and stability properties 
in this case are determined by the exchange process rather than the radiative cooling. 

Better fits are obtained by introducing another parameter describing a quadratic dependence on harmonic number, 
5j « (JqW^ -|- (5jQ (n — 1/2) -I- 5^ . This quadratic correction is very small in all examined cases, but its presence affects the values 
of and 5^ such that 5^ ^^lo sooner and more readily in the high-eg cases, as shown in Table ^. With this fit 5^ ~ I'^io 
by the point of eg = 100. Without the 5^ term, the sequence of modes approaches the behaviour of a doubly-open pipe much 
more slowly in Es- 

4.3 Stability properties 

The most unstable systems investigated are those with bremsstrahlung cooling only (ts = 0). Varying the values of as and 
-f/ioi has little effect in the small-ts extreme, because most of the cooling is due to bremsstrahlung radiation, which occurs in 
the dense region near the white dwarf surface, where the electron and ions have nearly reached equilibrium regardless of their 
initial two-temperature conditions close to the shock. For most realistic parameters with low or modest ts, the fundamental 
mode is stable. The fundamental is unstable only in cases of extremely low ijjci, modest or high (Js, and low ts. Otherwise the 
second harmonic is the lowest potentially unstable mode. 

Except when cyclotron cooling and two-temperature effects are both very strong (when -i/iei is very small and ts is large) , 
for given conditions of (o-s,i/'ci), increasing ts stabilizes every mode in a monotonic manner. Under these usual conditions a 
mode which is unstable for some value of ts is also unstable for all cases of lower es; and furthermore an increase in Es does 
not causes a decrease in a mode's stability. 

For a given change of ts the stability term 5-^ of each mode progresses in a way that is individual to the mode and 
apparently independent of the stability behaviour of the other modes of the sequence. Some modes experience more rapid 
reduction in 5-^ for a given increase in es, therefore they tend to stabilize earlier or quicker than other modes, if all else is 
equal. 

The transition of a particular mode from the unstable [5-^^ > 0) to stable (Jp, < 0) regions of the complex frequency 
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eigenplane occurs at some threshold value of ts that is individual to the mode and dependent on the (asjipai) conditions. The 
lowest modes usually stabilize at lower ts because they tend to be in or closer to the (5j^ < region even in the bremsstrahlung- 
cooling only case. However there often are exceptions, because of the the modes' rates of stabilization, dS^/dea, depend on 
Cs and are individual to each n. There are many cases of {as,ipBi, Es) where a mode n is unstable but a higher mode n + 1 is 
stable. Therefore for some cases of (cts, V'ei) there are modes which can never be the lowest unstable mode for any value of eg. 

It is instructive to examine collectively the stability values 5^^ of the sequence of modes at some set of parameters 
((Ta, i/Jci, fs). In the cases of lower eg there is a general pattern in which each successive harmonic is less stable than the lower 
neighbouring mode; 5^ increases with harmonic number n. The increases of between consecutive modes is smaller in the 
ranges of higher n. As eg increases this trend weakens, because some modes stabilize fast enough to overtake neighbouring 
lower modes. When eg is very high, the modal stability eigenvalues, Sj^, are scattered in a negative range of values, except 
for the lowest modes which tend to deviate further towards greater or lesser 5^. The introduction and strengthening of 
two-temperature effects enhances the departures from the general pattern. 

The eigenmodes for various eg under the one- and two-temperature conditions are compared in Figure |l|. Decreasing ipd 
from indefinitely high values (single-temperature fiow) to very low values (the two-temperature conditions strong) has little 
effect on the stability of the modes if bremsstrahlung cooling is the dominant process (tg is small) . In intermediate ranges of 
Eg, with cyclotron and bremsstrahlung cooling comparable in efficiency, strengthening the two-temperature conditions tends 
to stabilize most modes. When eg is very high, the infiuence of two-temperature effects on the stability of individual modes is 
pronounced. Introduction of two-temperature effect in high eg tends to destabilize most modes, and the pattern of variations 
of 5^ from one mode to the next becomes very different from the one-temperature case. For some {aa,ipei), two-temperature 
effects cause one or more of the lower modes to remain unstable even for very high eg, or even to destabilize a mode as eg 
increases past some threshold. 



4.4 Transverse perturbation 

The stability properties and frequency behaviour of the oscillatory modes were investigated in the presence of transverse 
perturbations, under the two-temperature conditions and with explicit competition between bremsstrahlung cooling and the 
second cooling (cyclotron cooling). The eigenvalues of the fundamental and first overtone were traced through a range of k 
for particular choices of ((jg, V'oi, eg). 

The special case k = is the absence of transverse perturbations in which the oscillatory behaviour is the same as that 
described in subsections above. In the limit of large k all of the modes become strongly stable, and the frequencies approach a 
sequence like a doubly-open pipe, oc 5ji~,n, Between these extremes, the efi^ect of transverse perturbations on the frequency 
and stability of the accretion column is specific to each mode. For a given mode with harmonic number n, the interesting 
behaviour takes place up to a threshold k « 2n-f 1. As k increases, there is a growing region of the eigenplane within which the 
modes n ^ (k — l)/2 show the simple behaviour like the large-K extreme. The modes outside the region experience individual 
and non-simple disturbances as k is varied, reverting to the standard a — behaviour in the large-n modes. 

For each mode, the instability term experiences a number of maxima as k varies. For a mode with harmonic number 
n, there are n of these maxima. Increasing eg afi'ects the values of these maxima in (k, 5^). For low eg the maxima are most 
distinct, but one or more maxima may become narrow and shallow almost to the point of vanishing as eg becomes very large. 
(See Figure 5, for example.) In the fundamental mode, the single maximum is at about k « 3/2 for low eg (see the eg = 
restricted case in Imamura et al. 1996) and gradually moves towards lower k as eg increases. In the first overtone, there is 
a maximum at k = and another one in the vicinity of k ~ 1 — 4. As eg increases, the former maximum lowers in and 
narrows, whilst the second maximum moves to lower k. When a maximum of instability (other than k — 0) migrates to smaller 
K with increasing eg, the tendency is for the peak to increase until some critical value of eg is attained, after which the 
peak decreases slowly with increasing eg. These thresholds depend upon n, ipd and o"s. The threshold is less than eg = 100 
when CTg is small. In instances of very high eg a mode may be stable for all k, except when (jg is large or i/^ci is small. In 
all cases studied, enormous values of eg J; 100 are required to stabilize modes other than the fundamental. For example, the 
fundamental is stabilized completely by eg = 100 in Figures ^, ^, but not in Figure ^ where ipd is too small. 

For all modes under all conditions, the oscillatory frequency 5j increases monotonically with k. In the presence of transverse 
oscillations, consecutive modes are no longer unique in 5-^ : it is possible for modes n and n -I- 1 to have the same but different 
transverse component k. However the modes do not become degenerate; there is no case in which different modes have the 
same k and S^. For each mode the rise of Sj with k is smooth and is a function of eg. Typical values of the gradient dS^/des, 
for given k, are greater for higher eg. For sufficiently high k and eg, 5^ reaches a plateau for each mode. (See the left side of 
Figures |-|). 

4.5 Summary 

We have found: 

(i) When eg is large, two-temperature effects are very important in determining the frequency and stability properties of 
the accretion fiow; when eg is small, the two-temperature conditions are ineffectual. Increasing eg decreases the frequency 
spacing of the modes. Varying eg alters the instability of each mode at a rate that depends on the harmonic number n, and 
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Figure 1. Complex 5 = <5p, + j<5j cigcnplancs, with contours of the logarithm of the perturbed velocity eigenfurictiori evaluated at 
the white dwarf surface. The modes are the points where the perturbed velocity goes to zero, marked with crosses. The left column 
represents the one-temperature cases; the right column shows two-temperature cases with 4'ci = 0.5 and CTs = 0.5. The top row shows 
bremsstrahlung-cooling only shocks, with eg = 0; in the middle row bremsstrahlung and cyclotron cooling are comparable €s = 1; and in 
the bottom row cyclotron cooling dominates Eb = 100. 
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a = 2.0, 13 = 3.85, a =0.50, f„ = 0.50 a = 2.0, /S = 3.85, a.= 0.50, f„=0.50 
1.01 I 0.051 




Figure 2. Oscillatory part and stability term of the eigenfrequency for the n = 1 mode, with varying transverse wavenumber k, in 
the presence of two-temperature effects. Squares, triangles and diamonds represent the bremsstrahlung plus cyclotron cooling cases 
Eg = 0, 1, 100 respectively. 



a = 2.0. (3 = 3.85, a =0.50, f„=0.50 a = 2.0, ,8 = 3.85, a.= 0.50, ^„=0.50 
2.0 I 0.10 1 




Figure 3. Same as Figure ti but for the n = 2 mode. 
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a = 2.0, /? = 3.85, fT,= 0.20, f„=0.50 a = Z.O., 13 = 3.85, a =0.20, Tp.=Q.5Q 

1.01 0.05 1 




Figure 4. Same as Figure Q but for the n = 1 mode, CTs = 0.2, i/iei = 0.5. 

a = 2.0. /S = 3.85, (7,= 0.20, •^,,= 0.50 a = 2.0, /S = 3.85, a.= 0.20, ^„=0.50 
2.0 I 0.10 1 




Figure 5. Same as Figure but for the n = 2 mode, as = 0.2, ipcl = 0.5. 
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a = 2.0, 13 = 3.85, a =0.50, 'i/'„ = 0.10 a = 2.0, ,3 = 3.85, ct.= 0.50, f,.= 0.10 




Figure 6. Same as Figure ^ but for the n = 1 mode, as = 0.5, 4'cl = 0.1. 

a = 2.0. /S = 3.85, a=0.50, f,,= 0.lO a = 2.0, /S = 3.85, a,= 0.50, f„=0,10 
2.0 1 0.101 




Figure 7. Same as Figure Q but for the n = 2 mode, CTs = 0.5, 4'ei = 0.1. 
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the two-temperature flow parameters {osjtpci)- Except when eg is very large and CTs is small, the 5j, of each mode only ever 
decreases with an increase in eg. 

(ii) The effect of varying the electron-ion coupling is insignificant in the bremsstrahlung-cooling dominated regime (es 
small). Decreasing ?/)ci, making the electron-ion exchange inefficient, tends to destabilize the accretion column in all modes, 
and increases the frequency spacing of the modes slightly. 

(iii) Varying the ratio of electron and ion pressures at the shock (jg affects the frequency spacing of the modes in non-simple 
ways. For some choices of the global parameters a given change in cts increases the frequency spacing; in other cases it causes 
a decrease. Decreasing as usually reduces the instability of modes, and this is accentuated when is high. 

(iv) Transverse perturbations yield 5^ maxima in k that can destabilize modes that are stable to purely longitudinal 
oscillations. If CTs is not too small then these maxima may enter the stable regime again for sufficiently great ts. 



5 DISCUSSION 

Most previous perturbative analyses of radiative shocks (e.g. Chevalier & Imamura 1982, Bertschinger 1986, Houck & Chevalier 
1992, Toth & Draine 1993, Dgani & Soker 1994) represented the cooling function as a single-power-law term, and set the 
indices to various values to mimic the presence of different physical cooling processes. Such formulations do not adequately 
describe systems where more than one cooling process is significant, especially where these processes differ greatly in their 
innate influences over thermal instability. The realistic interplay between the influences of a stabilizing cooling process and a 
destabilizing cooling process was explicitly included in Saxton et al. (1997) and Saxton et al. (1998). This condition is retained 
in the present extended treatment, and by taking the one-temperature limit we recover the earlier results. 

Our investigation considers the additional effects of unequal electron and ion temperatures, as in the work of Imamura et 
al. (1996), which considered accretion onto non-magnetic white dwarfs. Ours combines the two-temperature treatment with 
the effects of the joint presence of bremsstrahlung and cyclotron cooling. Therefore our analysis is applicable to radiative 
accretion shocks on white dwarfs in magnetic cataclysmic variables with radiative cooling timescales possibly comparable to 
or faster than the electron-ion energy exchange. 

The accretion column in a mCV is not geometrically extended in the vertical direction (see e.g. Cropper 1990), and 
the field-aligned flow suffers no transverse mass flux (unlike the case of Dgani & Soker 1994). Our the cases of nonzero 
correspond to magnetic fields much stronger than those investigated by Hujeirat & Papaloizou (1 



1998), and therefore the field 



can not develop large transverse kinks that would suppress oscillations by the mechanism of Toth & Draine ( 1993 ). The 
altitude of the shock above the white dwarf surface is negligible compared to the size of the white dwarf, so the variation of 
gravity in the accretion column, studied by Houck & Chevalier (1992), is unimportant. Therefore we do not investigate effects 
of non-planar geometry (as in Bertschinger 1986). Transverse perturbations, which are manifest as corrugation of the shock 
structure and oscillation, are studied, (as in Imamura et al. 1996). 

Except when the energy exchange of electrons and ions is very weak, increasing the efficiency of the cyclotron cooling Es 
(i.e. strengthening the ambient magnetic field) stabilizes each mode monotonically. The modes stabilize at independent rates 
in Es, meaning that some of those modes which are highly unstable in the weak-field regime may actually stabilize at low 
thresholds of Es. The detailed behaviour of the mode stabilization in eg depends on the parameters a a and ■i/'oi. 

As in Saxton et al. (1998) the instability of a particular mode does not imply the instability of all higher modes, when more 
than a single cooling process is significant. This characteristic of competition between cooling processes persists in the general 
two-temperature and non-longitudinal oscillation cases. If cyclotron cooling is significant relative to bremsstrahlung cooling, 
the electron- ion exchange process determines the sequence of 5-^, but the frequency spacing of the quantized 5^ sequence 
depends mainly on ta. Cases of bremsstrahlung-cooling dominated shocks are similar and their properties are scarcely affected 
by two-temperature effects. 

The basic sequence of mode stabilities in the bremsstrahlung-cooling dominated case has a stable fundamental and each 
higher mode less stable than its predecessor. Strengthening cyclotron cooling tends to stabilize all the modes, but does so in 
a fashion that is individual to each mode, so that deviations from the trend develop as eg increases (especially in the presence 
of two-temperature effects). For low eg these stability deviations are slight. For greater Es, the 5-^ sequence of the modes is a 
less simple pattern depending on the parameters of the shock pressures cts and electron- ion exchange ipei- 

Two-temperature effects dominate the properties of the modes when ts is very high, and beyond some threshold in eg, 
which depends on (cts, i/'oi), two-temperature effects have a destabilizing influence. If the electron- ion exchange is sufficiently 
weak, the exchange process dominates the oscillations and some modes that would be stable at lower eg become unstable 
again at higher eg. 

The non-simple order in which the modes are stabilized means that the thermal instability properties of accretion shocks 
in mCVs are not simply a trivial result of a relationship between oscillatory, cooling and energy exchange timescales. The 
competing effects must be considered in detail. However we note that our analysis proves modes to be unstable, but does not 
necessarily prove stability. Minorsky (1962, chapter 14) uses topological and perturbative analytic arguments to consider the 
general conditions of stability of non-linear oscillatory systems. This issue was also addressed by Lyapunov in his analysis of 
stability of dynamical systems (Lyapunov 1966, see section 1.5.3, Leipholz 1970 for another treatment of Lyapunov's theorems 
on stability). In summary, the linear terms dominate in the limit of small amplitudes, i.e. every mode that is linearly unstable 
remains unstable in a full non-linear treatment, although higher-order contributions may cause a stable limit cycle at larger 
amplitudes. 
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6 CONCLUSIONS 

We investigated plane-parallel radiative shocks in which bremsstrahlung and cyclotron cooling may be fast compared with 
electron-ion energy exchange. Parameters and boundary conditions are chosen to be appropriate for accretion onto the surface 
of a magnetic white dwarf. The cooling function was an explicit sum of contributions from thermal bremsstrahlung (which 
destabilizes the flow) and cyclotron cooling (which has a stabilizing tendency). 

The relative efficiency of the cyclotron emission es is varied. Stability of the shock to longitudinal and non-longitudinal 
perturbations was investigated. The ratio of electron and ion pressures at the shock as, and the rate of electron-ion energy 
exchange tpei were also varied. Variations of the exchange, pressure and cooling parameters, (us, V'ei, es), affect the stability and 
frequency of each mode in a manner that is individual to the mode and the conditions of the stream. Increasing eg increases 
the stability of each mode, until a threshold value in es is exceeded. Beyond this point, increasing es destabilizes a mode 
slowly. This threshold is very high except when ^pci is very small. 

The introduction of unequal electron and ion temperatures has little effect on the modes of a bremsstrahlung-cooling dom- 
inated accretion column. In these cases, the higher harmonics tend to be less stable to oscillations. When bremsstrahlung and 
cyclotron cooling are comparable in strength, the pattern of the modes' stabilities becomes less regular, and two-temperature 
effects tend to stabilize the shock. For a cyclotron-cooling dominated shock all of the modes are stable in the single-temperature 
limit, but introduction of two-tompcrature effects can make them less stable. In some cases of small tpci one or more of the 
lower modes may be unstable even for high eg. Unlike the bremsstrahlung emission, cyclotron cooling is strongest near the 
shock, where the difference of electron and ion temperatures is greatest. Therefore the slightest presence of two-temperature 
effects means that the stability of a cyclotron-cooling dominated shock is dominated by the energy exchange between the two 
fluid components rather than the cooling function. 

Many cases exist where the mode n is unstable while the mode n -|- 1 is stable, and some higher modes are unstable. This 
is a general characteristic of radiative shocks with competing stabilizing and destabilizing cooling processes, and is preserved 
in the generalization to conditions with two-temperatures and non-longitudinal perturbations. 

The oscillatory parts of the dimonsionless cigonfroquencies follow a sequence which is approximately regular and linear. 
In the singlo-tcmpcraturo case this sequence resembles the modes of a pipe open at one end: 5j ~ ^^^{n — 1/2) + S^^. When 
two-tcmpcraturc effects are present the limit of large es causes the frequencies to approach quantization like a doubly-open 
pipe 5; w Sj^n, because of the flxed velocity condition at the white dwarf surface becomes less signiflcant than the electron- ion 
thermal disparity throughout the column. The frequency spacing of the modes, S^^ , always decreases as es increases, but tends 
to increase when the electron-ion exchange is weaker. The frequency offset is small for low and it generally increases 
with increasing eg; in the limit of cyclotron cooling and two-temperature dominance the offset approaches —> ^Sjq, giving 
the doubly-open pipe behaviour. 

Introducing transverse perturbations can destabilize modes that are stable to purely longitudinal perturbations, except 
when cyclotron cooling dominates or two-temperature effects are strong. The modes become non-unique in S^, for different 
K. In the limit of large k all modes are stable with oscillatory parts quantized like a doubly-open pipe <5j oc n, regardless of 
the manner of their quantization under the same parameters in the absence of transverse perturbations. Each mode has a 
number of instability maxima in k, and these maxima evolve as the parameters {as,tpei, es) change. When the pressure ratio 
at the shock cxg is small, the maxima of 5^ are stabilized at lower eg. When the electron- ion exchange is weak (■i/'ei is small), 
the maxima may continue to destabilize at high ta- 
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